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Abstract. 

In this paper we study the evolution of a spherical matter overdensity in the context 
of the recently introduced Galileon field theory. Our analysis considers the complete 
covariant Lagrangian in four dimensions. This theory is composed by a potential 
and a standard kinetic term, a cubic kinetic term and two additional terms that 
include the coupling between the Galileon and the metric, to preserve the original 
properties of Galileons also in curved space-times. Here we extend previous studies, 
which considered both the quintessence and the cubic terms, by focussing on the role 
of the last two terms. The background evolution we consider is driven by a tracker 
solution. Studying scalar perturbations in the non-linear regime, we find constraints on 
the parameter of the model. We will show how the new terms contribute to the collapse 
phase and how they modify physical parameters, such as the linearized density contrast 
and the virial overdensity. The results show that the Galileon modifies substantially 
the dynamics of the collapse, thus making it possible to observationally constrain the 
parameters of this theory. 

Keywords: Modified gravity, Galileon, Spherical Collapse, Cosmological Perturbations 
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1. Introduction 

The discovery that the Universe underwent a phase of accelerated expansion at late 
times, through the study of the distance-redshift relation of type-la Supernovae (SNIa) 
[1, 2, 3], opened a new scenario for theoretical cosmology: the possibility to live in 
a Universe whose dynamics is presently driven by a component responsible for an 
"obscure" repulsive force, which has been dubbed Dark Energy (DE). Such a component 
should fill 74% of the energy budget of the universe, and it can be obtained by either 
just considering a non-zero cosmological constant term (ACDM model). This model fits 
very well observational data, but, up to now, it is impossible to give a physical meaning 
to the tiny value of A required to explain dark energy. Thus, cosmologists explored 
alternative theories by e.g, modifying the Einstein-Hilbert action: 



where M p \ represents the reduced Planck mass. Models that have been proposed are 
scalar-tensor theories [4], f(R) gravity (for a review see [5]), massive gravity (see [6]), 
Brane- World models (e.g. [7]) and others. 

Recently, a new class of theories was introduced by Nicolis et al. [8] , the so-called 
Galileon [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]. This model was 
constructed as an effective field theory, which is based upon and aims at extending the 
Dvali-Gabadadze-Porratti model (DGP) [26]. It is interesting because it solves ghost 
instabilities, which plague DGP, and has a screening mechanism that allows to satisfy 
the bounds coming from solar system experiments. To avoid the appearance of ghosts 
it is important to keep the equation of motion up to second-order in time-derivatives. 
Unfortunately, in the original model, this property was respected only in flat space- 
time. The works by Deffayet et al. [27, 28] found a way to generalize Galileons to 
curved space-time. To do this, it is necessary to add some extra terms which couple 
the scalar field with curvature terms. The result is a scalar-tensor theory in which the 
action, in flat space-time, is invariant under Galilean symmetry (<9 M — > <9 M + b^). 

Even though in this paper we study the effects of the late-time cosmic acceleration 
produced by the scalar field, adopting the spherical collapse model, it is worth 
mentioning that the importance of the Galileon field also relis on the fact that it can 
inspire some "inflationary-like" model, e.g. [29] (even though this theory does not 
respect the galilean symmetry). 

In this paper we will use the background evolution given by the tracker solution 
found in [30], which ensures a de Sitter (dS) stable point. While it is shown that the 
cosmology of the Galileon lets the universe expands accelerating at late-times, in this 
paper we will show that at short distances we can satisfy solar system constraints via 
the Vainshtein mechanism [31, 32], see also [33] for a discussion in the most general 
second order scalar-tensor theory. This mechanism can work in massive gravity, but 
also in different contexts. An example is DGP theory, which possesses a Vainshtein 
radius defined by r v = [r s r1) 1 ^ (where r s is the Schwarzschild radius of the source, 




(1) 
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r c is a coupling constant which defines the crossover scale between a 5-dimensional 
Minkowsky space and the embedded 4-dimensional space-time). Even if the Galileon 
does not already have a well-defined Vainshtein radius, we will show how to recover a 
valid definition of it. In fact, as in DGP, instead of a massive graviton, this mechanism 
can also work using non-linear self-interaction terms of the scalar field (as D0(V0) 2 ). 

The spherical collapse model (e.g. [34, 35, 36]) studies the evolution of a spherical 
Dark Matter (DM) overdensity to explain the formation of cosmic structures. We 
will use the top-hat approximation, taking into account the energy non-conservation 
problem noted in [35] . This problem affects theories with a time-dependent dark energy 
component, and it can substantially modify the virialisation process. 

The paper is organized as follows. In Section 2 we define the action we are assuming 
and we obtain the equations of motion. In Section 3 we briefly review the background 
evolution of a Friedmann-Lemaitre-Robertson- Walker universe (FRLW) following a 
tracker solution found in [30]. In Section 4 we study scalar perturbations, both in 
the linear and non-linear regime. We also study the Vainshtein mechanism, and discuss 
the existence of a solution for the Galileon field in the non-linear regime. In Section 5 we 
study the dynamics of a spherical top-hat matter perturbation. In Section 6 we discuss 
our main results. In Appendix A and Appendix B we give some useful functions. 

Throughout the paper we adopt units c = K = G = 1, except where explicitly 
indicated; our signature is (— , +, +, +). 



2. Action and Field equations 



Let us start with the covariant action for the Galileon model non-minimally coupled to 
the metric [30]: 



/ 



Mix 

-f R+ 2 



1 5 



+ / d x C M , 



(2) 



where q are dimensionless constants. We consider Cm as the Lagrangian of a pressurless 
perfect fluid with density p. The five Lagrangian densities for the scalar field are: 

d = M 3 (f) (3) 

A = (V0) 2 (4) 

£ 3 = (D0)(V0) 2 /M 3 (5) 

£ 4 = (V0) 2 [2(D0) 2 - 2<t>,^ - i?(V0) 2 /2] /M 6 (6) 



, 5 = (V0) 2 [(D0) 3 - 3(D0) 4>,^ + 20/0/0/ + 



60 ;M 0^0^G//M 9 



•mi- > (7) 
where M is a constant with dimensions of mass, and we defined its value as M 3 = 
M p iH% s . iJjs is the value of the Hubble parameter H(t) in a FRLW universe at the 
de Sitter fixed point. Indeed, as we will see, [30] found a tracker solution that ends 
at a stable point called "de Sitter point", at which the energy density of the scalar 
field dominates. C\ can be understood as a potential term and for this reason we set 



Spherical Collapse in covariant Galileon theory 



4 



Ci = 0, since we are interested in analyzing the contribution of the new kinetic terms 
(the case in which a standard minimally coupled scalar field is introduced in the field 
equations was already studied in [37]). Moreover, with this choice we can employ the 
tracker solution given in [30], that is not admitted if c\ ^ 0. £2 is the standard kinetic 
term. £ 3 comes directly from the decoupling limit of DGP theory. £ 4 and £ 5 provide 
the full generalization of an action containing at most second derivatives with respect to 
Galilean shift symmetry in a flat space-time. The coupling between and the curvature 
tensors are required to construct a Lagrangian free of third or higher-order derivatives 
in the equations of motion. 

Varying this action with respect to the metric g^ v and the scalar field we obtain 
the equations of motion. For the metric: 



G 



T 



(m) 



+ T 



where 



(8) 



T 



(<t>) 



(0 



(9) 



the terms being listed in Appendix A. Instead, varying with respect to the scalar 

field, we obtain 



£q£« = 0, 



i=i 



where £W are also listed in Appendix A. 



(10) 



3. Background evolution 



From Eqs. (8) and (10) we can study the background evolution in an expanding FLRW 
universe with scale factor a(t). Calling = <p(t) and p = p m (t) + p r (t), the background 
scalar field and background matter and radiation density respectively, the field equations 
read 



and 



where 



3M 2 { H 2 

3M 2 U 2 + 2M 2 } H = /', 



P<P + Pm + Pr , 

I' 



p r /3, 



C 2 



3H<p + (j) 



+2H(f) + 3H(f) 



3c 3 

M 3 
15c 5 

~~ M 9 



3H 2 (f) + H<p + 2H<p 



M 6 



3H 2 (f)+ 



2i3 



3H 2 (j) + 3H(j) + AH(f)\ = , 



2 r 



3c 3 

M 3 " 



45c 4 



H 2 ^ + , 



21c, 



315 



2M 6 

3C4 



2 

-^ 9 H^[5H4> + 2(H 2 + m, 



M 3 



2M 6 



M 9 

3 [8#0 + {3H 2 + 2H)4>} + 



(11) 
(12) 

(13) 



(14) 



(15) 
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are scalar field density and pressure, respectively. 

As in [30], to study the background we work with the new variables 

ri = <j> dS H dS /(j>H) , r 2 = (0/0 dS ) 4 /n , fi r = p r /(3M 2 tf 2 ) , (16) 

where </>ds is the time derivative of the scalar field at the dS point. At this point Eqs. 
(11) and (12) becomes: 

c 2 x 2 dS = 6 + 9a - 12/3 , (17) 
c 3 x 3 dS = 2 + 9a - 9/3 , (18) 

where a^s = 0ds/(-^ds^pi)- These equations give two conditions for the coefficients c 2 
and c 3 . We also set a = c^x\ s and /3 = c 5 x dS ; therefore our free parameters become 
a, (3 and Xds- For simplicity, the assumption Xds = 1 will be often used in the rest of 
the paper. An approximation we have done is H dS ~ H , where H is the value of the 
Hubble parameter today. 

As we already mentioned, [30] found a stable tracker solution {r\ = 1), which drives 
the universe expansion from the radiation-dominated epoch (r2 < 1, fl r = 1), through 
the matter-dominated epoch (r 2 = 1, Q r 1), until the dS point (r 2 — 1, Q r — 0). 
Note that along r\ — 1, = p^j {?>M 2 X H 2 ) = r 2 . Following this solution, Eqs. (12) and 
(10) with our new variables can be written as 



r 9 = 



2r 2 (3 - 3r 2 + tt r ) 



7r 2 ) 



(19) 



l + r 2 ' r l+r 2 

where primes denote differentiation w.r.t. = In a. In Fig. 1 we show the numerical 
solution of these equations with boundary conditions fi ro = 4.8 • 10~ 5 and f2 Ao — 0.74, 
where Q ro and fl\ are the density parameter values today, for the radiation and the 
dark energy component, respectively. These equations cannot be solved analytically; 
however we have found two analytic functions that approximate the numerical results 
with an accuracy better than 1.2% at redshift z < 21: 



r 2 (A0 ~ 1 + 



[i-n Ao ) 2 l-n 



Ao 



2Q 



Ao 



/ 4e 6/V + 



(1 - ^A )^ 



Ao 



Ao 



• e 



-6N 



(20) 



and: 



n r (N) ~ 2^"" (i - ^A + V / 4^A e 6iV + (l-fiA )^ • (21) 

To study the stability of the solution ri(N) = 1, Eqs. (11), (12) and (13) can be 
expanded at linear order in perturbations 5r\, 5r 2 and 5Q r . Thus, it can be obtained: 

9 + Q r (N) + 3r 2 (N) 



5n'(N) 



2(l + r 2 (iV)) 



-5n(N), 



which reads: 



Sn(N) = 5ri(0)exp 



^ 9 + ^(^0 + 3^(^0 
u ' l 2(l + r 2 (iV0) 



(22) 



(23) 



/o is a finite integration constant, and this relation proves that any solution that 
approaches r±(N) = 1, finally reaches it. Indeed, in the rest of the paper, we shall 
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Figure 1. In the figure we show the evolution of tt r (red line) and (green line), 
functions of N = In a. 



suppose that at least after the matter-dominated epoch the evolution of the universe 
can be described by 8r\ <C 1. 

In [30], the authors also find constraints on the parameters a and (3 (assuming 
^ds = !)• These constraints follow from the requirement of ghost avoidance. They 
study scalar (S) and tensor (T) perturbations, expanding the action Eq. (2) at second- 
order in perturbation theory (see [38, 39], for the complete procedure), finding conditions 
for the sign of the kinetic term (Q$ and Qt) and the squared sound speed (c s 2 and c T 2 ). 
Thus, in every epoch we have four conditions that must be satisfied. Reminding that 
a and (3 are constants, we can find a region of parameter space where no ghost modes 
exist. This area is bounded by the analytic functions 

' a > 2/3 
a < 2/3 + 2/3 
a < 12^-9/3-2 
a > 12/13/3 + 10/13. 



(24) 



4. Cosmological perturbations 

In this section we study the evolution of scalar perturbations on sub-horizon scales. 
Our work focuses on the dynamics of a spherically symmetric perturbed metric. Let us 
choose the conformal Newtonian gauge, 

ds 2 = -(1 + 2^)dt 2 + a 2 (t)(l + 2$)5 ij dx i dx j . (25) 

Perturbations of the energy density and the scalar field are given by 

p(x,t) = p (t) + 8p(x,t) (f)(x,t) =(f>o(t) + <P&t). (26) 

In the following we will drop the suffix "0". In this regime there are two 
valid approximations that simplify the field equations. The first one is the sub- 
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horizon approximation (9(V 2 $/a 2 ) ^> 0(H 2 &). The second one is the quasi-static 
approximation, which allows us to neglect time derivatives of perturbations compared 
with space derivatives, assuming we are working with non-relativistic matter at short 
distances. 



4-1. Linear perturbation theory 

Replacing physical gradients with comoving gradients, at linear order Eqs. (8) and (10) 
become (V denotes a spatial gradient): 

(2M 2 + 2 7 i(f)) V 2 $ = -bp + 72(f) V 2 y? , (27) 
(2M 2 + 373(f)) V 2 $ + (2M 2 + 2 7 i(t)) V 2 ^ = 374(f) VV , (28) 



and 



75(f) VV + 72(f) V 2 ^ + 374(f) V 2 $ = , (29) 



where 7»(f) are functions of the background, whose explicit form is given in Appendix 
B. 

It is important to note that one of the differences between these equations and those 
for the kinetic braiding model studied in [36] is the presence of an anisotropic stress in 
the RHS of Eq. (28). 

Manipulating Eqs. (27) and (28), we obtain the modified Poisson equation 

2 



(2M 2 +0 2 7i 
2M 2 + 3 73 



-V 2 ^ = 5p- 



+ ^71 

72 d75 2M 2 +3 7 3 



VV (30) 



Using Eqs. (29), (27) and (30), the differential equation for the evolution of the 
scalar field takes the form 

V 2 <p = A(t)6p(t,r), (31) 

where 

A { t) s 2 7 2 (t)7 7 (t)- 37 4 (f)7a(t) 

7 2 (f) 2 7r(f) - 7 5 (f)7e(f) 2 - fry^fr^freW 

with 7 6 (f) = ^2M 2 j + 2 7i(f)) and 77(f) = (2M 2 t + 373(f)). Considering a spherically 
symmetric object of radius Rs, we can easily integrate Eq. (31) to obtain an analytic 
expression for the evolution of the scalar field. Defining m(t,r) = 4rr J Q r dr'r /2 5p, we 
obtain 

£ = «t )mi t,r) + C 

dr Airr z r z 

where C is an integration constant that, outside the source, can be viewed as an increase 
in M s = m(t, Rs)- While this term is present in tp', it does not enter in V 2 <p, so that the 
gravitational potential is not affected by our choice of C. Therefore, for our purposes 
we can set C — 0. 
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4-2. Vainshtein mechanism in the linear regime 



The Vainshtein mechanism works by screening the effects of the scalar field on the 
gravitational potential at small distances, so that one can satisfy the constraints coming 
from solar-system tests, while preserving the accelerated expansion of the universe on 
cosmological scales. The difference between this mechanism and the Chamaleon one is 
that the first also works by using non-linearities of the perturbations to this aim. At 
large distances (r ^> ry, where ry is the Vainshtein radius of the source) linear terms of 
the scalar field become dominant, while for r <C ry non-linear terms become dominant 
(these terms will be shown in Eqs. (40), (41) and (42)). This is called "self-screening 
effect". A discussion about the magnitude of the Vainshtein radius (r v ) of a spherically 
symmetric source will be given later (Sec. 4.4). 

A first approach is to study within the linear approximation the contribution of 
the scalar field to the gravitational potential. Recalling Eq. (30), to have a qualitative 
knowledge that outside the Vainshtein radius the scalar field drives the late time cosmic 
acceleration, we have to compare the contribution of the gravitational with the scalar 
field intensity [40]. Indeed, our request is that the two are comparable: 



ip'(r) 



1. 



(34) 



#'(r) 

It can be shown that the above ratio is a monotone function, which starts from ~ 
during the radiation-matter-dominated epoch. At the dS point, recalling Eq. (31) with 
Xds = 1, we obtain 

(35) 



Taking into account the region in the plane (xds = 1,13, a) bounded by the no-ghost 
condition (24), it can be shown that the magnitude of the last ratio at the dS point is 
bounded by 







A(t dS ) 




1 


*'(r) 


dS 


47T 




24ttM p1 (2/3 - a) 



W2n 



< 



ip'(r) 



#'(r) 



< +oo 



(36) 



dS 



This result means that the contribution of the scalar field at the dS point on scales 
r ^> ry is always important, and the importance can be set choosing proper values for 
a and f3. In particular we can find a couple (a, (3) which satisfies Eq. (34). 

With Eq. (31) we can write the modified Poisson equation (30) in a more convenient 
form: 



where: 



V 2 # = 4vrGV5p(t,f), 



75(077(0 + 9 7 4 (tf 



4tt [6 72 (t)74(07 6 (0 - 72(0 7t(0 + 75(076(0 ] 



(37) 



(38) 
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The modified gravitational constant assumes the value of the Newtonian one during the 
radiation-matter-dominated era, while it is 

G 



= W^W) (39) 

at the dS point (when xas = 1)- The limit Xds gives us the usual GR result 
G^t^s) — G. Instead, the limit Xds — > oo gives G^t^s) 0, which means, as expected, 
that the effective gravitational constant becomes small w.r.t. the Newtonian one 
(G oc MpJ 2 ). The plots in Figs. 2, 3 and 4 show that we can vary the asymptotic value 
of G^ as we desire, to obtain, in principle, any reasonable model for the late time cosmic 
acceleration. The difference between the three graphs is the value of the parameter 
Xds, which sets the contribution of the Galileon field at the dS point. This result also 
agrees with the expectations of Eq. (36), quantifying the effective contribution of the 
scalar field at large distances on observables quantities. Of course, these results do not 
represent any realistic model, we are only interested here in investigating the range of 
possibilities offered by the Galileon theory. Moreover, astrophysical and cosmological 
constraints on the Galileon model have just started being considered [41, 42, 43, 44]. 




N 



Figure 2. This plot shows the evolution of G^, with Xds = 1, in different 
cases. The values for (a, ft) are: (—1,-0.55), blue dashed line; (—0.45,-0.4), red line; 
(—0.2,-0.2), green line; (—0.55,-0.4), blue solid line; (0.1,-0.1), red dashed line. 



4-3. Non-linear evolution 

When perturbations grow, Eqs. (27), (28) and (29) must be replaced by fully non-linear 
ones. Neglecting time-derivatives of perturbations and assuming that the characteristic 
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Figure 3. The same as in Fig. 2, but with x^s = 0.3. 




Figure 4. The same as in Fig. 2, but with x^s = 1.2. 



scale of the perturbation is well within the Hubble radius, we obtain 

(2MJ + 2 7 i(£)) V 2 $ = -dp + 72 (t)VV + li(t) [(VV) 2 + 
-V>V>] + m (t) [(VV) 3 + 2V>V>V fe ^+ 



3 ■ 



3VVV>V J '^] - -<p\(t) [vVv 2 $ - v^v^] 
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(2MJ + 373(0) V 2 $ + (2M 2 + 2 7 i(O) V 2 ^ = 3 74 (t)V 2 ^ + 



+ 3 V2 (t) [(W) 2 - v>v>] - l<j> 2 m(t) [VW 2 * - V^V>] . (41) 

Eq. (29), instead, takes the form 

75 (t)V 2 v9 + 72 (t)V 2 ^ + 3 74 (t)V 2 $ + 773(0 [(V 2 ^) 2 - V>W 

- 774(0 [(V 2 ^) 3 + 2V>V>V fc ^ - 3VVV>V>] + 

+ 2 7l (0 [V 2 y?V 2 ^ - V l ^V>] + 6?7 2 (0 [V 2 v?V 2 $ - v^v>] - 

- i^V(0 [v 2 ^v 2 $ - v^w^] + 3?7i(0 [(vV) 2 v 2 ^- 

-2VVV>V j ^ - V 2 W>V> + 2V>V\y?V fc ^] = , (42) 

where the r]i(t) functions are listed in Appendix B. 

Eqs. (40), (41) and (42) are more complicated than in the linear case, however, 
assuming spherical symmetry, they are in fact integrable. The boundary values of the 
perturbations can be determined by resorting to the physical meaning to these fields. 
For example, from GR we know that the physical solution of the Poisson equation is 

Gm(t, r) 



*Gfl'(*,r) 



(43) 



Recalling the definition of the mass function, m(t,r) = 4rr J Q r dr'r' 2 5p(t,r), if there 
are no singularities at r = for the density perturbation, this relation tells us that 
^Gfi'(^O) = (to violate this limit we have to choose 5p(r) oc r~ n , with n > 3). 
At small scales we want to recover GR, so the physical meaning of \l/(t,r) should be 
that of gravitational potential (\l/'(t,r) ~ ty GR '(t, r)). Indeed, the natural assignment is 
^'(t,0) = 0. The same argument applies to &(t,r ->■ 0) ~ -^ GR '(t,r ->■ 0). Instead, 
the scalar field and its perturbations are not directly observable quantities, so we have 
to choose the correct boundary value by mathematical arguments or by its effect on 
measurable physical quantities. Like in Eq. (33), at r — > there should be some 
divergent term. However, the same reasoning used in the linear case allows us to consider 
ip'(r — > 0) finite. 

Integrating Eqs. (40), (41) and (42) for a spherically symmetric object, we obtain 



76" 



r 



m(t, r) 
4 7rr 3 



72" 



2 7 i— + 2771^ - 2771 — 30 771 — ^— 44 



77 V 7e— = 374 h 677; 



12 



30 ^1—5- 



(45) 



75 h 72 h 374 h 27/; 



12 



r r 

\3 



r 



277, 



/3 



+ 



+ 67/; 



+ 1277: 



30V 



= 0, (46) 
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Note that we have not yet analyzed the case in which the scalar field has a boundary 
value ip'(r = 0) finite, but different from zero; to do this we have to impose a physical 
condition. From Eq. (45), we can write 



-74 



<f'(r) = _j ± _ fjh _ ^(r) Sgn( 74 ) 
r Ar] 2 Ar]2 r Ar] 2 

+ (2M 2 + 3 73 ) ^ + ^2 (2M p 2 , + 2 7 i) 



+ 



1/2 



3 v pi r 3 v - 7 r 

here we have chosen the solution which matches the linear one (33) when r — > oo. 
Without any loss of generality, the metric perturbations can be written as 

V(r) = *' G «(r) [1 + MO] = [1 + MO] 

$'(0 = * G «(0 [1 + MO] = [1 + MO] . 



(47) 



(48) 
(49) 



In this case, ^Gi?( r ) can be understood as the gravitational potential generated by a 
perturbation in the ACDM model. When r <C ry, 5^ and 5$ have to be small by solar- 
system constraints ((5^,(5$ < 10~ 3 ), so we can treat them as small perturbations. In 
this limit, at first order, Eq. (47) becomes 



74 gVj^Oj Sgn( 74 )/(t,r) 
4r]2 4?7 2 ^ 4^2 



Ar] 2 r 



+ 



Sgn( 74 )fr' GR (r) 

12r] 2 f(t,r)r 
Sgn( 74 ) 77 ^ GR (r) 
/(t,r)r 



30 V 2 ^^-30 2 74^i + 4^2 



M^O , 



-M*,r) + 
M*,0 + 
(50) 



where 



/(«, r) 



873^2 



;2 *Gfl(0 



8 - 2 *' Gfl (r) 
3 r 



1/2 



(51) 



From Eq. (50), we are now ready to choose a reasonable boundary value for (p'(r). 
It is sufficient to suppose that neither 5^ nor (5$ diverge in the limit r — > 0, to show that 
(f/{r ->■ 0) = 0. 



Vainshtein radius 

Having obtained the non-linear equations of motion, we are now ready to investigate 
the radius at which non-linearities become important. The simplest way to estimate ry 
is to plug-in the linear solutions into the non-linear equations, and estimate when the 
non-linear terms become comparable with the linear ones. First, considering Eq. (44), 
from the quadratic term we obtain 
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To solve this equation, we need to know the matter density profile. However, using Eq 
(33), in the general case we find 

li(t)A(t) 



[m(t, r) + Am(t, r, ryj] , 



(53) 



27T 72 (t) 

where Am(t,r, ryj = 4-7T J^ Vl dr'r' 2 5p. The interior solution for a top-hat profile leads 
to an r-invariant equation. The simple consideration is that, depending on the epoch 
and on the choice of the background parameters, we can have this region all inside or 
all outside the Vainshtein region. Instead, outside a source of mass M s we find (defining 
Ry = r v (R)) 



R 



Vi 



R 



271 

72 



A(t)M s 



AnR 3 



The same procedure for the cubic term leads to 

3 



R 



v 2 



R 



2m 

72 



A(t)M s 



AttR 3 



(54) 



(55) 



Comparing the two Vainshtein radii we see that they are comparable. This means that 
we have an exterior linear region, but, when we enter the non-linear one, quadratic and 
cubic terms can both dominate. Indeed, the contribution derived from the terms C4 and 
c 5 influences in a non-negligible way the scalar field profile. This also proves that at 
sufficiently large distances we recover the predictions of the linear theory, discussed in 
Sec. 4.1. 

Other three important Vainshtein radii, coming from Eqs. (45) and (46), are 



(56) 



R v ^ 


\ 3 = 2m 


A(t)M s 


R , 


1 74 


4ttR 3 


Rva^ 


\ 3 = 2r?3 


A(t)M, 


R , 


' 75 


AttR 3 



and: 



Ry* 
R 



2ru 

75 



A(t)M s 
AttR 3 



(57) 



(58) 



Here we have neglected non-linear interactions which couple tp with $ and 
because they produce results analogous to the previous ones. The Vainshtein radius 
can be set as Ry = Max(R Vi ), where % — 1, .., 11. It is straightforward to prove that 
R v (t ->■ -00) ->■ +00, while R v (t ->■ +00) = f(a, f3, x dS )M s /(AicM p iH^ s ), where / is a 
generic function of the background parameters. This result agrees with the predictions 
of [35] and [36]. 



4-5. Galileon field evolution 

In this section we study the Galileon field evolution, starting from Eqs. (44), (45) 
and (46). These are three algebraic equations in ^'(r), $'(r) and y?'(r), so it is 
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straightforward to obtain a sixth-order polynomial equation in <p'(r) (to simplify the 
problem we will work under the assumption that x ds = 1): 

m' 6 a/ 5 (z/ 4 u?' 3 

^ + \l(t)^r + Ht)^ + (WWm + A 4 (t)) ^ + (X 5 (t)S m + 
12 I 

+A 6 (t)) ^ + (A 7 (t)5 m + A 8 (t)) ^ + \ 9 (t)5 m + A 10 (t)5 m 2 = , (59) 

where Aj are background functions, combinations of 7, and 77,. From Eq. (59) it follows 
that <//(r) has six branches of solutions. What is the correct one? Remembering the 
Vainshtein effect, we want that the physical solution reduces to Eq. (33) at large 
distances. Of course, this condition cannot be verified analytically, but it is sufficient to 
choose between the real solutions of Eq. (59). 

Are we sure that, for a given couple (a,/3), Eq. (59) has at least a couple of 
solutions during the whole evolution of the universe? Obviously this condition is not 
sufficient to ensure the existence of the physical solution, but it is a necessary condition. 
In the linear regime the existence of a physical solution was proved in Sec. 4.1, thus the 
problems can be inside the Vainshtein radius. As proved in Sec. 4.4, at small distances 
non-linear terms become dominant for the evolution of the scalar field. In particular, 
instead of Eq. (59), we can work with the equation 

m' 6 (n' 5 U ,A 

+ Ai(f)^- + A 2 (*)^- + \ w (t)5 m 2 = 0. (60) 

Also in this case we do not have an analytic solution for the scalar field; however Eq. 
(60) gives new constraints on the allowed region in the parameter space (a, (5). 
Consider a function like 

f(x)=x e + Ax 5 + Bx 4 + C, (61) 

where A, B, C 7^ are real coefficients. The RHS of this equation has the same form 
as Eq. (60), after the substitution ip'(t,r)/r — > x. It was demonstrated that there is no 
analytic method to find a solution for f(x) = 0, when f(x) is a fifth or higher degree 
polynomial. However, since 

lim f(x) = +00 , (62) 

it is sufficient to require that a minimum of this function is < 0, to be sure to have at 
least a couple of real solutions. The points which satisfy f'(x) = are 



5 / 25 2 

*iA3 = x^ = --A±^—A*--B. (63) 

The zeros of Eq. (59) can be understood as six perturbative terms around x,. Let us 
assume that, for the purpose of this section, these perturbations are small. The set of 
parameters for which f(x) =0 has at least a couple of solutions, which are given by 

/(xi)<0 V /(x 4 )<0 V /(x 5 )<0. (64) 
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Substituting our background functions into the parameters A, B, C, we must pay 
attention to the dependence on t, because the previous inequalities have to be hold true 
Wt. The first one becomes 

4 



H ds M P1 ) 



4 + 



H 12 M 6 
1440 4 /3 2 

' H ds Mpi 

It can be proved that this condition is verified Vt, when 

a < 4/5 

a < 5.22/3 + 1.93 
a < -3.73/3 + 4.83. 



5« + 42/3(^) 



8 m {tf < 0.(65) 



(66) 



These relations were obtained evaluating the above expression at some critical times, 
when f(xi) results maximized/minimized. We were able to do this because /(0) takes 
a simple form, but this is not the case for f{x$) and f{x§). In fact, the form of these 
functions at the points x^ 5 is 

f(x 4>5 ) = c-— (±5A + 725^-96^) 4 (5A 2 - 245+ 

±AV25A 2 



96B 



(67) 



In our case, the parameter C depends on the matter-density perturbation, so the 
inequalities which follow from the above expression have to be evaluated in two distinct 
cases. The first one is when the density term dominates on the other terms (the analysis 
is the same as in /(0) < case), the second one when it is subdominant. The latter 
case involves more complicated expressions for the parameters a and (3, so we were only 
able to solve it numerically. Combining these results with the no-ghost condition given 
in [30], the constraints on the parameters a and (3 become 

a > 2/3 
a < 2/3 + 2/3 
a < 4/5 
^ a < 5.7/3 + 2.62, 

and are represented in Fig. 5. 



(68) 



5. Spherical Collapse 

In this section we will restrict our analysis to a top-hat matter configuration 

*>-{r* :ii^H^ (rlRf (69) 

The mass 5M is the total mass of the density perturbation Sp, while M = 
4/3 n (po + Sp) R 3 . The two masses are related by 

5M = M , (70) 

1 + 6 
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a = 4/5 L. 




a = 2/3 + 2/3 


■ / 
~-/\ 


-0.8 -0.6 -0/ -0.2 


\ 0.2 0.4 

^ \ 


/ -0.5 


: \ 

a = 2 p 


/ / \ _1 ° 




a= 5.1/5 + 2.6 




-1.5 





Figure 5. In f/iis figure we show the allowed region in the plane (/3, a) obtained by 
mixing the no-ghost conditions and the conditions for the existence of the scalar field 
in the non-linear regime. 

where S = Sp/po is the density contrast. 

To study the dynamics of a spherical matter perturbation we need the well known 
equation 

a k*2 

which follows from the non-linear continuity and the Euler equation for a pressureless 
fluid of non-relativistic matter in a top-hat configuration [34]. Eqs. (40), (41) and (42) 
tell us that, inside a top-hat density perturbation, ^/'(r) oc r, which means that V 2 $ 
will be r-in dependent. Indeed, a top-hat profile, remains a top-hat profile during its 
whole evolution despite the non-validity of Birkhoff 's theorem. 

To solve Eq. (71), we have followed [34]; here we briefly summarize the main steps. 
Assuming the total mass conservation, R 3 po (1 + 5) = const., Eq. (71) can be rewritten 
in terms of R 

— = H 2 + H — -V 2 ^. (72) 
R 3 

From this equation we can distinguish all the sources that affect the collapse dynamics: 
H 2 + H contains the contribution of the background (matter and dark energy), while 
V 2 \I/ contains the contribution of matter and scalar field perturbations. Using N = In a 
as a time variable and defining 
R a 
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where Ri and a« are the initial radius of the perturbation and the initial scale factor, 
Eq. (72) becomes 

y" + I wy , -( 1 + I p)y = -l(y + eN ~ N ') v2 * > ( 74 ) 

where a prime denotes differentiation w.r.t. N. The density contrast is 

8 = (l + 6 i )-(e N '- N y + iy 3 -l. (75) 

Eq. (74) can be solved numerically setting the initial conditions. From Eq. (73) we 
know that yi — and y[ = — 5^/(3(1 + <5j)). Supposing that the perturbations start 
growing linearly during matter-dominance, the linearization of Eq. (71) can be solved 
analytically. The growing mode is 5 oc a, so 5' = 5, thus the second initial condition 
becomes y[ = —Si/3. We also set a* = 10 -5 , while the initial density perturbation is set 
to collapse exactly at clq — 1. 



5.1. Virialisation 

The Virial Theorem states that a stable system must satisfy the relation 

W + 2T = , (76) 

where 

T=\j SxpJ = ^Mi? 2 (77) 
is the kinetic energy (the last equality holds true for a top-hat profile), while 

W = _ fd 3 xp m (x)x -W = -^V /V-r 3 ^i (78) 

y iv 3 ~* Jo dr 

is the trace of the potential energy tensor. As in the previous equation the last equality 
holds true only for a top-hat profile. \&i(r) denotes each component that contributes to 
the total gravitational potential. 

Usually energy conservation is used, but, as noted in [35], for a time-dependent 
dark energy model, energy is not strictly conserved. So, during the collapse phase, the 
virial radius can be estimated as the radius at which the virial condition (76) is satisfied. 

Important quantities that can be extrapolated from the dynamics of the collapse 
are the linearized density contrast 8 C , and the virial overdensity: 

Kir = = [1 + S(R vir )} ■ f^vA 3 . (79) 

Pcollapse V ^vir J 

5.2. Numerical Results 

5.2.1. Case (3 = 0, Xds = 1- This is the case in which the fifth term of Eq. (2) gives 
no contribution. Eqs. (44), (45) and (46) become simpler. In particular, the modified 
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Poisson equation reads 



V 2 ^ = 3^ m # d 2 s a- 3 a; 4 2x " - 5 + 



3H 2 s x 2 [2x 4 {2 + a) 



a 



-2 + 15a)] - 36a(2x 4 + 3a)H if' 



H 2 s M pl (2x 4 + 3aY 



+ 



\2ax 2 {2x 4 - 3a) <^ n 
~ H 2 s M 2 j2x 4 + 3a) 2 ' V 2 
with x = H/Hds, and <// is a solution of: 



with: 



if' 3 ip' 2 ip' 
a.i ■ — + a 2 ■ — =- + (0:3 + «4<5) • h a 5 5 = , 



«i = 4ax 2 (4x 8 + 24x 4 a - 45a 2 ) 

a 2 = 2M p i [F 2 s x 2 (4a; 4 (2 + 3a) (x 4 + 6a) - 9a 2 (2 - 21a)) + 

+6a (4x s - 24ax 4 - 45a 2 ) H 

a 3 = - 2^ S M 2 ! [i^| s x 2 [2x 8 (2 + a) + x 4 (-4 + 8a + 21a 2 ) + 
+a (2 - 21a + 45a 2 )] - [4x 8 (2 + 3a) + 27a 2 (-2 + 5a) + 

+12a; 4 a(-2 + 9a)] H 

a 4 = - 8e- 3n ^ s M 2 ,Q m x 4 (2x 4 - 3a) a 

a 5 = - e- 3n 4M p 3 ,fi m x 2 [F 2 s a; 2 (2x 4 (2 + a) + a(-2 + 15a)) + 
-12a (2x 4 + 3a) H . 



(80) 

(81) 

(82) 
(83) 



(84) 
(85) 

(86) 

Of course, among the solutions we want the one that reduces to Eq. (33) when 
r > r v . 

Although this is a particular case, it is interesting to show the role of £4 in Eq. (2). 
In Fig. 6 we have plotted the solution of Eq. (74) for various a. It should be noted that 
modifications w.r.t. the ACDM model are present during the collapse phase. This is, as 
expected, an effect of the increasing contribution from the scalar field. In Tab. (1) we 
show the values assumed by the linearized density contrast and the virial overdensity. 



5.2.2. Case a = 0, x^s — 1- In this paragraph we analyze another particular case, the 
one which shows the role of £5, Eq. (2), in the dynamics of the collapse. Compared 
to the previous paragraph, when (5 7^ Eq. (59) cannot have an analytic solution. By 
the parameter conditions, Eqs. (24) and (68), — 1/3 < (3 < 0, so, to investigate the 
parameter region in which (3 > we need to set a > 0. 

The dynamics of the collapse is shown in Fig. 7, while the linearized density contrast 
and the virial overdensity for various (3 can be found in Table (2). It is important to note 
that the onset of the fifth term in Eq. (2) plays a crucial role in the virialisation process. 
In fact we can see that varying the parameter j3 there is a substantial modification of 
A vir with respect to the ACDM model. 
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41 


914 
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a = 1/10 


2.243 


1.723 


0.537 


28380 


44 


0.899 


14500 


328 


a = 1/5 


2.272 


1.757 


0.527 


27930 


46 


0.884 


14850 


305 


a = 1/3 


2.300 


1.801 


0.515 


27470 


48 


0.863 


15430 


272 


a = 1/2 


2.327 


1.847 


0.504 


27020 


51 


0.836 


16150 


238 


a = 2/3 


2.345 


1.882 


0.495 


26680 


53 


0.812 


16710 


215 



Table 1. Here we show numerical results of physical interesting quantities in the case 
j3 = 0, Xds = 1 for various a 




Figure 6. In the figure we plot the solution of Eq. (74), in terms of the normalized 
radius R/Ri of the top-hat perturbation, when (5 — and Xds = 1- The initial density 
for each model is shown in Tab. (1). 

5.2.3. Case a / 0, ^ 0, xas = 1- This is the most general case, despite the 
assumption Xds = 1. Here we can evaluate the sum of the contribution of the terms £4 
and £5, Eqs. (6) and (7), in the whole parameter region defined by Eq. (68). It can be 
noted that as a and (3 grow we obtain a larger 6 C , thus it should be easy to remove a 
large piece of parameter space from the allowed region. 

6. Conclusions 

In this paper we have first reviewed the background evolution of the Galileon model, 
following the tracker solution of [30] . We have found two analytic functions that describe 
the evolution of the components of the universe at late-times. The peculiarity of this 
tracker solution is that it ensures a dS stable point independent of the parameters 
of Eq. (2). This assumption simplifies our equations, but it should also be easy to 



Spherical Collapse in covariant Galileon theory 



20 




Figure 7. In the figure we plot the solution of Eq. (74), in terms of the normalized 
radius R/Ri of the top-hat perturbation, when a = and x<is = 1. The initial density 
for each model is shown in Tab. (2). 
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46 
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= -0.1 


2.296 


1.780 
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27620 


47 


0.866 


15060 
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/3 = -0.2 


2.356 


1.857 


0.501 


26740 


52 


0.813 


16440 


225 


13 = -0.3 


2.412 


1.928 


0.484 


25980 


57 


0.769 


17150 


198 



Table 2. Here we show numerical results of physically interesting quantities, in the 
case a = 0, Xds = 1 for various j3 

generalize our work to a more general background evolution. Once c\ is set to zero, in 
Eq. (2) should remain only kinetic terms for the scalar field, thus the Galileon cannot 
be considered as a deviation from the ACDM model. It should work as a substitute of 
the cosmo logical constant, mimicking the effects of A on cosmological scales. 
Then we have shown that, in the linear approximation the scalr perturbations of a 
FRLW universe lead to a time-dependent gravitational constant G^it), that modifies 
the gravitational potential generated by a distant or, equivalently, small source. The 
results we give do not represent a realistic model, i.e. they are not required to satisfy the 
observational bounds, rather they are chosen on order to display what one can generally 
expect from this theory. 
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R/Ri 




Figure 8. In the figure we plot the solution of Eq. (74), in terms of the normalized 
radius R/Ri of the top-hat perturbation, when a ^ 0, (3 and x^s = 1- The initial 
density for each model is shown in Tab. (3). The values for (a,f3) are: KCDM blue 
dashed line; (—0.45, —0.4) red line; (—0.2, —0.2) green line; (—0.55, —0.4) blue thick 
line; (0.1,-0.1) red dashed line. 
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Table 3. Here we show numerical results of physically interesting quantities in the 
case a = 0, Xds = 1 for various (3 

The Galileon model can be successful because it possesses a Vainshtein mechanism, by 
which we can consider two distinct regions; the first one at large scales, where the linear 
approximation applies and the Galileon drives the cosmic acceleration, the second one 
where non-linearities are dominant. We have also shown how to recover a Vainshtein 
radius in agreement with the one of DGP and other simpler models. 
Even though the study of the perturbations in a highly non-linear regime can notbe 
completely analytic, we found some constraints, whose fulfillment allows Eq. (60) to 
have at least a couple of real solutions. 

The last part of this paper was devoted to the study of the collapse of a spherical top- 
hat matter perturbation. We have shown that the new terms £4 and £5 affect in a 
non-negligible way the dynamics of the collapse and the value of 5 C and A V i r . To study 
the virialisation process we paid attention to the energy non-conservation problem, 
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Appendix A. Components of the field equations (8-10) 



The terms of the stress-energy tensor of the scalar field read 

r (1 > =£mV* 



T 



(2) 



fiu 



T 



(3) 



T 



(5) 



-0;^0;z,+ (W) 1 



T (4) = 

/IV 



- X« [0;m 0;^ D - 0;{m 0;<4« ; ° + ; ° 0;«/3 ;/? ] 
1 



M 3 

2 
M6 



(A.l) 
(A.2) 
(A.3) 



:R<P;^;u (V0) 2 + 20 ; ^0^0 ; ^^ + 



-20 ;M „ ;a ;a/3 ;/3 + 20 ;{M ;/3 <f)^ a - ;M ; „ 0^+ 

+ Ra{» <P;v } <?* (V0) 2 - ^(V0) 4 + ;/la 0/ (V0) 2 + 

-9^ R a p ;a (V0) 2 + ;a ;/3 (V0) 2 + 

-2g^ ;a ;/3 ;/?7 W - ; ^ D0 (V0) 2 - 20 ;{M ;£/}a ;a D0+ 

+0 ;M ; „ (D0) 2 - ^ ;a/3 ;a/3 (V0) 2 + 2^„ ^ a ;a/3 ;/3 D0+ 



(D0) 2 (V0) ; 



(A.4) 
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[30: 
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7 A;" 



+ ^ 0^ 0^ (V0) 2 - -i^ ; , }/3 ^ (V0) 2 + 



2 

+^0 ;{M ;1 , }a r(V0) 2 + ^G 



4 
3 



(V0) 2 + 



--Rf ;Ut ; , }a ;/ * (V0) 2 - 30 ;Ata 0^ 0' a ;7 0^+ 

-30 ;MQ 0^ 0;^ (V0) 2 + 30 ; ^ ;a 0^ ;7/3 <P + 
-30 ;{M 0;, }a 0^ ;7/3 0^ + ^ 0; M 0;, 0^ (V0) 2 + 

+ ;M ; „ 0./ ;/37 0^ + ^ ;/i 0;„ D0 (V0) 2 + 

+ ;a /3 ;Q/? (V0) 2 + ^0 ;{M ; „ }q r 0;/3 7 0^ + 

0;, 0;a/3 0^ 00 + |i2o7^ 4>. v " ^ ^ (V0) 2 + 
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— Rate 0;,} 0;/3 0^ (V0) 2 - \r^ 0;,} ^ 0^ (V0) 2 + 

~R ;M ; , D0 (V0) 2 + ;Q D0 (V0) 2 + 

+3^ <P;a ;/3 0^ ;7 r 0^ + 3^ i? 7/3 ;Q ;7 0^ (V0) 2 + 
0; 7 0^ <^ 7 (V0) 2 - ;a ;7 0^ (V0) 2 + 

+ ^ i?a 7 /3r 0^ 0^ (V0) 2 + 9{lv 0^ ;a7 ;/?7 (V0) 2 + 

+ 30;^ 0;^ 0^ D0 - 30 ;/1 „ <?« ;a/3 0^ D0 + 

+30 ;{M ; , }Q ;/3a ;/3 D0 + 30 ;Ma 0/ D0 (V0) 2 + 

-\%» Rap ; ° ;/? D0 (V0) 2 + ;a ;/3 D0 (V0) 2 + 

(Q0) 2 (V0) 2 - 3 -<j>., { , <f>.„ }a ;Q (D0) 2 + ^ ; , (n0) 3 + 

-\g,u 0^ ;a/3 0^ 0; 7 r 0^ " ^ 0;a/3 0^ U(f> (V0) 2 + 

-3^ ;a ;7 ;7/3 ;a/3 D0 + ^ ;a ;a/3 ;/3 (D0) 2 + 

+^(D0) 3 (V0) 2 
The terms appearing in the equation of motion for the scalar field read 

P = ^ (A.6) 

e (2) = - D0 (A.7) 

e (3) = ^ [-(°^) 2 + ^ + ^ r v \ (A.8) 

e (4) = ^ [2R r ; ^ ; " - 8iU ; ^ ; ^ 0./ - 2R, V 0^ ( V0) 2 + 

-4^ 0^ 0^ 0^ - 40./ 0./ 0^ Q + i2 (V0) 2 D0+ 
+4i^ ; ^ ;i/ D0 + 60 ;/u , ; ^ D0 - 2(D0) 3 ] (A.9) 
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+6i? Q/3 ^ V ;M Q ; / + 0^ 0./ (V0) 2 + 
+60 ; ^ ; „ a ;Q/3 ; ^ - 80'^ <j>, va 0./ D0+ 
+12^ 0^0^ ;MQ 0^^ - 6i2^ ; ^ (V0) 2 D0+ 
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-3(0;^ + <f>" r ;7 ° </>/ - QR^aufS ^ 0^ □</>+ 

+12R^ <jP r ; / ^ + 3iW r v ^ (V0) 2 ] (A.10) 
Appendix B. Background functions 

Background functions involved in the linear perturbation theory: 
7i(t) = 3 (a - 2x dS /3) 

H dS M p\ 

2 



72 (t) = (2 + 9a - 9/3 - 12x dS a + 15x 2 dS /3) 



73(0 = - QEj4 ^2 « + ^ 



3^ d 4 s M 2 V ^ 2 s M pl 



74W = - nTj2 u ( 2a ~ 3x asf3) XdS + 772— 



7 5 (t) = - 6 - 9a + 12/3 - 26x dS a + 4x dS (2 + 9a - 9/3) + 2Ax 3 dS /3 + 
+ 2 [2 + 9a - 9/3 - 6x dS (a - x dS /3)] -A- 

Background functions involved in the non-linear dynamics: 

2/5 12 



^(*) = orrl ^2 I " - 6/3- 



3^ d 4 s M 2 I ^ 2 s M pl 



2 + 9a - 9/3 - 6 (a - x dS /3) ( x dS + 



^S^PI 



2 1 

= ^ 4 ^2 « - 2/3- 
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